Statistics of Intense Turbulent Vorticity Events 
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We investigate statistical properties of vorticity fluctuations in fully developed turbulence, which 
are known to exhibit a strong intermittent behavior. Taking as the starting point the Navier-Stokes 
equations with a random force term correlated at large scales, we obtain in the high Reynolds 
number regime a closed analytical expression for the probability distribution function of an arbitrary 
component of the vorticity field. The central idea underlying the analysis consists in the phase-space 
restriction to a particular sector where the rate of strain and the rotation tensors can be locally 
regarded as slow and fast degrees of freedom, respectively. This prescription is implemented along 
the Martin-Siggia-Rose functional framework, whereby instantons and perturbations around them 
are taken into account within a steepest-descent approach. 

PACS numbers: 47.27.Gs, ll.15.Kc 



The intermittent nature of turbulent flows, long ago 
discovered by Batchelor and Townsend l], has been an 
issue of main concern in fluid dynamic research. De- 
spite the great investigation efforts taken so far, boosted 
in the last two decades by considerable experimental 
and numerical progress, there is not yet a clear theo- 
retical understanding of intermittency (see Ref. 2] for 
an overview). In fact, the sudden and intense vorticity 
bursts observed at high Reynolds numbers can be ad- 
dressed in general along two distinct - but not necessar- 
ily incompatible - viewpoints. Strong fluctuations may 
arise as a consequence of local vorticity amplification by 
stretching, induced on its turn by strain defined at larger 
scales. Alternatively, high peaks of enstrophy and dis- 
sipation are also detected whenever a vortex coherent 
structure is advected across the measurement position. 
Such an interpretation ambiguity was emphasized in the 
recent experimental work by Zeff et al. 3], where some 
support was given for the dominance of the former, local 
mechanism of vorticity bursting. 

Intermittency is frequently characterized in homoge- 
nous isotropic turbulence by the statistical properties of 
the velocity derivatives sij = diVj, which, as expected on 
general grounds, have their second and higher moments 
diverging as some power of the inverse viscosity A 
relevant problem is to find the profile of the probability 
distribution function tails (pdf tails) of Sij. Exact an- 
swers have been found in the simpler realm of the Burg- 
ers model, where the presence of shock waves leads to 
a bifractal form of intermittency j^. It turns out that 
in general the velocity derivative pdf tail converges in 
a pointwise sense to C|s|~^/^ 0, Q, for s < 0, in the 
vanishing viscosity limit v (for s > the pdf de- 
cay is faster than gaussian). In contrast, much less is 
known in the case of three-dimensional incompressible 
turbulence. Actually, just a few attempts to describe the 
viscous asymptotic expressions for the pdfs of the veloc- 
ity derivatives are available so far, devised on t he g rounds 
of phenomcnological multifractal ideas 0, 0, Ii(tI |. which 



stress the dissipative role played by small scale vortex 
structures. 

The main result of this letter, to be obtained with the 
help of functional techniques, is an analytical expression 
for the tails of the vorticity pdf p{uj) in the high Reynolds 
number limit, where lu is an arbitrary cartesian compo- 
nent of the vorticity vector uji — eijkdjVk- Our starting 
point is the effective modeling of turbulence, as described 
by the stochastic Navier-Stokes equations, 

dtVa + VfsdpVa = -9qP + vd'^Va + fa , daVa = . (1) 

The gaussian random force fa{x, t) has a vanishing mean 
value and the two-point correlator Dap{x — x')5(t — t'), 
with 



Dap[x - x') = 5al3DoG-iq^{~\x - x'\^ / L^) 



(2) 



Assuming the Kolmogorov cascade picture energy is 
injected into the fluid at rate Z?o, within large length 
scales of order L, and dissipated at the microscopic scale 
77 ^ v^l^ — > 0, where viscous effects come into play. 

Let r = dxaVa be the circulation evaluated at time 
t = 0, for a circular contour c of radius i? — > 0, counter- 
clockwise oriented and centered at the origin in the xy 
plane. The vorticity pdf p{w) may be written in terms of 
the characteristic functional Z{X) = (exp(zAr)) as 



lim — 

fl— 27r 



dX e^p{-iXujR'^)Z{X) 



(3) 



The radius R enters in the formalism as a regularizing 
parameter to handle some singular expressions which ap- 
pear in the course of computations (i? ^ is taken at 
fixed v). The Martin-Siggia-Rose functional approach 
jllj l may be now recalled to put Z{X) into a p ath inte- 
gration form. Up to a normalization factor |l2l |. we have 



Z{X) = J DvDvDPDQexp{iS) 



(4) 
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where the Martin-Siggia-Rose action is given by 



S= Sx dtC + XT 



(5) 



related to the non-local lagrangian density 

C = VaidfVa + VfjdfSVa - l^d'^Va + daP) + QdaVa 



+ I d x'Dapix - x')v0{x',t) 



(6) 



Of course, the limit T — > oo is intended in Q. How- 
ever, we will resort in practice to a temporal domain 
of size T = kL^/j/, having the order of the viscous de- 
cay time of smooth large scale fluctuations (introduced 
in the forthcoming instanton computations), where n is 
some Reynolds number dependent dimensionless param- 
eter. Free boundary conditions are allowed at t = — T, 0. 

As an attempt to single out the phase-space sector 
where strong fluctuations of T develop, we perform the 
analytical mapping A — > —iX in (jSJ, taking afterwards 
the large A limit. The saddle-point technique may then 
be employed to set non-perturbative field configurations 
(instantons) which would solve the variational equation 
SS = 0, and hopefully yield the dominant asymptotic 
behavior of Z{X), directly related to the form of the 
vorticity pdf tails. These computational steps, which 
can be carried out for the analysis of pdf tails of gen- 
eral observables, constitute the essence of the instanton 
method successfully applied to problems like the 

random advection of a passive scalar and Burgers turbu- 
lence 0, Q, 0, 01 . Unfortunately, the implementation 
of the instanton approach in the three-dimensional in- 
compressible situation is plagued with non-perturbative 
difficulties, as previously noted in Ref. [l^. It is some- 
what simple to understand why the saddle-point action 
that would be obtained along the above lines cannot lead 
to the correct vorticity pdf tails (or to the pdf tails of 
any linear combination of velocity derivatives). Since the 
saddle-point equations are invariant under the transfor- 
mations 

i,'^ = hv^ , P' = hP , Q' = h^'^Q , A' = /lA , (7) 

where h is an arbitrary scaling parameter, it follows that 
the saddle-point action has necessarily the general form 
S = A^/^/(A~^/^i^), yielding, in the high Reynolds num- 
ber limit, ln[p(w)] —oj^ at the rigth tail (that is, uj > 0), 
which is a very unlikely expression under the light of nu- 
merical and experimental findings 0, 0, 0, 1201] . 

One might conjecture that the functional approach 
would work if fluctuations around the instanton were 
taken into account, possibly to all orders in the pertur- 
bation series, to overcome the X^^^ dependence in the 
saddle-point action. However, this is an extremely hard 
task to accomplish. A more computationally oriented 



point of view is to hypothesize that a family of instan- 
ton solutions could be defined from the start, so that 
the overall contribution of the associated saddle-point 
actions would lead ultimately to physically meaningful 
results. Indeed, we have found that the latter option is 
the most promising one. 

The hint on how to set up a family of instanton config- 
urations comes from the hypothesis that the small scale, 
intense bursts assumed to dominate the inviscid vortic- 
ity pdf tail are driven by the rate of strain provided by 
larger eddies, which according to standard ideas, fluc- 
tuate on slower time scales [2lj. In fact, several tur- 
bulence models depict vorticity enhancement along the 
energy cascade as a local process which takes place in a 
static background defined by an external straining flow 
[2^. Among the manifold scenarios proposed in the lit- 
erature, it is worth mentioning Lundgren's modeling of 
stretched vortex tubes with internal spiral structure [23j , 
intimately connected with the observations reported in a 
remarkable recent experiment by Cuypers et al. |24l |. 

The above phenomenological considerations motivate 
us to revisit the functional strategy initially founded on 
Eqs. H3I6(I . The point is to separate, in the integration 
measure, the slow degrees of freedom (the rate of the 
strain tensor) from the fast ones (the rotation tensor). To 
achieve this goal, we insert the identity 1 = / Da'^S[a^p — 
{daV/B + dfjVa)/2] in the integrand of Defining 



Z([cr''],A) = / DvDvDPDQDQ exp{iS) 



(8) 



a strain-dependent characteristic functional, associated 
to the modified Martin-Siggia-Rose action 

S = S+^J d'ldtQ^p [d^vp + dpv^ - 2<^) , (9) 

we write, after exchanging the order of integrations, the 
alternative (and exact) expression for the vorticity pdf in 
the conditional form, 

p{uj) = lim — P>(^" dAexp(-awE2)Z([CT"],A) . 

R^n 271 J 

(10) 

It is clear, thus, that instanton configurations may be 
naturally labelled (in a functional sense) by the rate of 
the strain tensor cr^^, if they are now found as solutions 
of the saddle-point equations derived from SS = 0, viz., 

A+Va - VpdaVf} + dp{vpVa + Qap) + iX- = , 

OVct 

A^Va + VpdpVa + daP + iDap W/J = . (11) 

Above, we have A± = 9* ± vd'^, Q*^ = Qap + QSap, and 
= e30o.^6irj_ - R)6{z)dit) , (12) 
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where r± = ^Jx^^^\^]p- . As can be easily checked, the 
saddle-point Eqs. Hll() are invariant under (|7I), replacing 
Q by Q*^ and adding (ct^^)' — h^^'^cr^fj to the set of 
transformations. At high Reynolds numbers, the modi- 
fied saddle-point action will have then the functional form 
5 = A3/2/[A-i/V]. 

It would not be an easy matter to solve exactly Eqs. 
However, as a general and valuable feature of the 
instanton method, it is not necessary to know the de- 
tailed global form of solutions of the saddle-point equa- 
tions 0, Q, 0, 01 . In our particular problem, the 
small radius R of the integration contour for F, which 
is in the viscous range, indicates that a linearization of 
t) around x = would suffice. Also, as a first order 
approximation, the rate of the strain will be considered 
time- independent, as suggested from its role as a slowly 
fiuctuacting variable in the vorticity intensification pro- 
cess. 

While we are not worried with the global spatial be- 
havior of the velocity field, the instantons are required 
to have finite action, and to satisfy the boundary con- 
dition 0+) = 0, assuring that solutions for Va(x,t) 
will not blow up for t — > oo, once the parabolic opera- 
tor A+ is non-negatively eigenvalued. We may replace, 
equivalently, the Va{x,Q'^) = condition by 



(13) 



as it follows from the evolution equation for Va {x, t) 
and the former boundary condition itself. Taking now 
dpQap = VjsidaVp - dpva) in (inj, we get, for t<Q, 



A+Wq - VpdfSVa + VfjdpVa = , 



(14) 



which is a Helmholtz-like equation. It is in fact always 
possible to obtain a symmetric tensor solution for Q*^, 
closely related to the result that the support of Va{x, 0~), 
a singular ring, is backwardly advected in time by the re- 
versed velocity field, as implied by Eq. H14|l . It is not dif- 
ficult to prove that the backward advection of the singu- 
lar ring will generate a contour with bounded length (and 
hence a non-vanishing saddle-point contribution) for all 
t < if and only if erf 3 = cr|3 = 0, and the two eigenvalues 
of in the xy plane are positive. Therefore, it is worth 
noting that the functional integration in (|10|) turns to be 
effectively replaced by ordinary integrations over the lin- 
early independent matrix components afi, erf 2 and a^^ 
of the rate of the strain. Observe that (t|2 is not in the 
set of integration variables, since cr|2 + '^11 + '''33 = 0; 
according to the incompressibility condition. 

We focus the analysis of the saddle-point equations 
on the class of axisymmetric instantons, having in mind 
the axial symmetry of Z([a''],A). The general form of 
Z([cr*],A) will be found through a symmetrization pro- 
cedure, taking into account that Z([cr''],A = 0) is invari- 
ant under transformations of the full rotation group, and 



any eventual A-dependent correction factor will be invari- 
ant only under rotations around the z axis. In concrete 
terms, we are just taking advantage of the fact that Eqs. 
(|ll|l become fully axisymmetric and consequently more 
manageable if we write the rate of the strain tensor - 
the static background in the saddle-point equations - as 
~ 2(5q3(5/33), where a is an arbitrary posi- 

eti 

saddle-point configuration 



-'a/3 



five parameter, and (5^^ = e^a-y^siij- We get, then, the 



Q{x, i) = , P{x, t) = (cj2(t)/4 - a^)r\/2 - 2o?z'^ , 
Wq(x, t) = o^^f^xfi - Lo{i)ezai3Xjil2 , 

X f-i 

Va{x,t) = iXt3,f3a — 5{ri_ — i?exp(ai))(5(z) , (15) 



where is the Heaviside function, and 
7rDnAi?2 



uj{t) 



lL2 



■ exp(— 2a|i|) 



(16) 



is the instanton vorticity. The imposition of Wq(0, i) = 
in p5(l is guaranteed from the invariance of 5* under 
generalized galilean transformations , while the pres- 
sure P{x,t) is fixed up to some unimportant additive 
time-dependent function, associated to the large scale 
boundary conditions. In order to consider the effects of 
smooth fiuctuations around the instanton solutions, we 
take the fields given in H15|l and perform the substitu- 
tions ^ Va + 5va, Va ^ + Sv^, ctc. in the path 
integration ||SJ), keeping in S terms defined up to the sec- 
ond order in the perturbations. The integrations over SQ 
and SQ enforce the constraints daSvfs + dpSva = and 
daSva = 0, respectively. Therefore, still working in the 
linear approximation for the velocity field, we will have 
Sva{x,t) = ^apxp, with $a_3 = eajS-fip-y- Defining, fur- 
thermore, Vap{t) — J cPxXaSvp (an antisymmetric ten- 
sor, as it follows from the integration over 5P), we obtain 
the phase-space reduced expression 



Z([a^],A) = exp(- 



4 \ 2 



tt'DqR'^X 
2L'^a 



DvapD^ap exp(-Af ) , 
(17) 



which depends on the effective action 



1 



(18) 



The subsequent integration over Vaf) gives 



Z{[a%\) = exp(-^-^^^) / D4>a exp(-i/) , (19) 



where we have now 
L2 



H 



2Do 



di(0Q 



^[0^(O) + 0^(-T)] 



(20) 
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The functional determinant factor in H19(l can be ex- 
actly computed, from its relation to the problem of a 
set of three one-dimensional euclidean harmonic oscilla- 
tors [2^ . The term proportional to e — *■ 0+ in H provides 
a zero-mode regularization of the functional integration. 
We find, thus, taking together the saddle-point action 
and the fluctuations around the instantons, 



Z(H,A) = 



cxp(7r^i:'ofl4AV£V|3) 
\/det[cosh(crT)] 



(21) 



The rate of the strain tensor which appears in H21|) can 
be diagonalized through a simple rotation by an angle 9 
around the z axis. Let afi = p — q, cr|2 = P + 9 ^n.^ 
CT^g — — 2p be the matrix elements of in its diagonal 
form, where p and q are real parameters. The jacobian 
of the mapping {p,q,0) (cfi, (7^2, crlg) is 4|g|, and we 
write, on the basis of H10|l . 



i?' 

p(ijj) = lim — / dp I qdq 



de 



d\eyi^{-iXujR^)Z{[a%\) 



(22) 



In the large T limit the hyperbolic cosine in (jnil is well 
approximated by a simple exponential and the vorticity 
pdf may be written in a closed analytical form. We get 
the t-distribution 



p{uj) 



16(tj2 + c,(t^2^); 



(23) 



where g (x k and (w^) is the actual expectation value of 
the squared vorticity, to be taken from experiments or 
numerical simulations (we have ^ diouP' p{lo) = g{uj'^)/A). 
We expect Eq. (|23|l to apply only at the pdf tails and 
their leading order corrections, where \uj\ ^ y/ (w^). 

An excellent agreement holding around three decades 
of the vorticity pdf is found between the tails predicted 
by (|23|l . with g ~ 5.5 and (w^) ~ 10, and the results of 
numerical simulations at moderately high Reynolds num- 
ber [23, as shown in Fig.l. The relevant vorticity range 
is localized considerably beyond the usual gaussian core. 
It is an interesting problem, left for further investigation, 
to determine the Reynolds number dependence, if any, of 
the fitting parameter g. 

To summarize, we studied the form of vorticity pdf 
tails at high Reynolds numbers by means of the instan- 
ton technique. We took into account, as an essential phe- 
nomcnological motivation, the idea that the intermittent 
vorticity bursts observed in turbulent flows arise from 
the coupling of vorticity and relatively slow fluctuations 
of the rate of the strain tensor. A supporting comparison 
was established between the vorticity pdf tails given by 
^ and previous numerical results. 
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FIG. 1: The vorticity pdf tails taken from (solid 
lines) are compared to the results of numerical simula- 
tions by Farge et al. at the Taylor-scale Reynolds 
number R\ = 150. 



